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SUMMARY 

Lockheed  Missiles  and  Space  Company,  Incorporated,  has  performed  a  program 
of  spectral  imaging  observations  of  the  solar  atmosphere  using  the  Sacramento 
Peak  Vacuum  Tower  Telescope.  The  observations  were  obtained  with  Lockheed 
instruments  including:  an  active  tilt  mirror  for  image  motion  compensation; 
polarization  analyzer;  narrowband  tunable  birefringent  filter;  photoelectric 
CID  array  camera;  digital  video  image  processor;  and  a  microcomputer  controller 
Five  observing  runs  were  made,  three  of  them  with  the  entire  system  in  opera¬ 
tion.  The  images  obtained  were  processed  to  measure  magnetic  and  velocity 
fields  in  the  solar  photosphere  with  very  high  spatial  resolution  -  0.5 
arcseconds  on  the  best  frames.  Sets  of  these  images  have  been  studied  to 
address  three  scientific  problems:  (1)  the  relationship  among  small  magnetic 
flux  tubes,  downdrafts  and  granulation;  (2)  the  puzzling  flux  changes  in 
isolated  magnetic  features  in  a  decaying  active  region;  (3)  the  temporal  power 
spectrum  of  the  magnetogram  signal  in  isolated  flux  tubes.  Examples  of  the 
narrowband  images  are  included  in  the  report,  along  with  abstracts  and 
manuscripts  of  papers  resulting  from  this  research. 
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BACKGROUND 

Lockheed  Missiles  and  Space  Company,  Incorporated,  has  carried  out  a  program 
of  spectral  imaging  observations  of  the  solar  atmosphere,  using  the  Sacramento 
Peak  Vacuum  Tower  Telescope.  All  components  for  a  sensitive  photoelectric 
spectral  imaging  system  were  assembled,  utilizing  previous  developments 
sponsored  by  Air  Force,  NASA,  and  Lockheed  Independent  Research  funds.  These 
instruments  were  used  to  study  convective  flows,  hydrodynamic  wave  modes  and 
magnetic  flux  tubes  on  the  sun  with  high  accuracy  and  resolution.  Previous 
observations  by  Lockheed  personnel  had  concentrated  on  studying  the  three- 
dimensional  structure  of  magnetic  fields  in  the  photosphere  and  chromosphere 
and  on  learning  the  most  informative  diagnostic  techniques  which  are  possible 
with  tunable,  narrow-band  filter  observations.  This  work  had  revealed  a  rich 
variety  of  phenomena  which  could  be  studied  by  spectral  imaging.  However,  the 
time  dependence  and  dynamics  had  not  been  addressed  in  any  detail  because  of 
the  lack  of  an  automatic  filter  control  system  and  the  data  handling  problem 
with  large-volume  film  and  video  recording  of  raw  data.  Over  the  past  five 
years,  we  acquired  the  components  of  an  automated  digital  system,  largely  as 
engineering  models  for  flight  hardware  of  the  Spacelab  2  Solar  Optical 
Universal  Polarimeter.  Specifically,  we  operated  an  image  motion  compensation 
system,  polarization  analyzer,  tunable  filter,  photoelectric  CID  array  camera, 
digital  video  image  processor  and  microcomputer  controller.  The  video 
processor  derived  images  of  brightness,  velocity  and  magnetic  fields  from  raw 
CID  frames  and  made  video  movies  of  these  processed  images.  These  data  were 
analyzed  using  the  image  processing  systems  at  Lockheed  and  at  Sacramento 
Peak. 

INSTRUMENTATION 

The  major  components  of  the  spectral  imaging  system  are  briefly  described 
in  this  section.  Conceptually,  these  comprise  an  engineering  model  of  the 
Lockheed  Solar  Optical  Universal  Polarimeter  (SOUP),  which  will  fly  on  the 
shuttle  flight  Spacelab  2  in  the  summer  of  1985.  The  SOUP  instrument  is  a 
much  more  powerful  observing  tool,  and  the  engineering  models  have  only  once 
played  together  to  simulate  Its  operation,  in  October,  1984.  Nevertheless, 
the  engineering  models  represent  a  substantial  investment  of  effort  and  funds, 
and  they  have  obtained  unique  and  exciting  solar  observations  to  date. 

Subsets  of  them  will  continue  to  observe  the  sun  periodically. 

Tunable  Filter:  This  is  a  Lyot-type  birefringent  filter  with  a  75  mA  band¬ 
pass  at  5250  A.  It  can  be  tuned  to  any  wavelength  between  5000  and  6600  A, 
provided  an  8  A  blocker  is  available.  A  HeNe  laser  provides  wavelength 
calibration,  accurate  to  a  few  mA.  Once  solar  line  positions  are  stored 
initially,  they  can  be  recalled  with  arbitrary  offsets  from  line  center, 
again,  accurate  to  a  few  mA.  Optically,  the  filter  is  of  good  quality,  as 
evidenced  by  the  very  sharp  images  (0.4  -  0.6  arcsecond  resolution)  obtained 
with  it  since  1979. 

Pref liters:  A  wheel  containing  six  waveplates  provides  polarization  analysis 
in  right  or  left  circular  and  four  linear  polarizations.  An  identical  wheel 
containing  dielectric  interference  filters  blocks  the  spectrum,  passing  one  of 
eight  regions  centered  on  5173,  5250,  5576,  5876,  5896,  6302,  6328  or  6563. 
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Image  Motion  Compensation:  An  agile  tilt  mirror  located  at  a  telescope  pupil 
stabilizes  the  image  despite  atmospheric  wander  and  telescope  shake.  It  has 
piezoelectric  transducers  for  speed,  a  simple  integrator  controller,  and  a 
quad  cell  centering  an  image  of  a  sunspot  to  provide  the  error  signal.  When 
operated  properly,  the  system  stabilizes  the  image  of  the  sunspot  (or  the  mean 
position  of  a  spot  much  larger  than  an  isoplanatlc  patch)  to  a  small  fraction 
of  the  seeing  disk. 

CID  Cbmera :  This  SOUP  engineering  model  uses  a  GE  CID-11B  array  image 
sensor.  It  has  248  x  244  pixels  of  46  x  35  microns  each.  It  can  be  read  out 
non-destructively  in  0.5  seconds,  digitized  to  12  bits.  The  best  sensors  can 
attain  slgnal-to-noise  ratios  of  500-1000  and  are  linear  to  a  fraction  of  a 
percent  (after  a  weak  pre-exposure  flash  of  light).  Fixed-pattern  gain 
variations  show  a  checkerboard  pattern  at  the  5%  level.  In  nearly-collimated 
monochromatic  light,  large,  smoothly-varying  fringes  are  seen  at  the  20% 
level.  Both  of  these  patterns  cancel  out  to  a  fraction  of  a  percent  in 
magnetograms,  which  are  ratios  of  an  image  difference  to  an  image  sum.  Thus, 
real-time,  relative  differential  photometry  of  high  sensitivity  can  be  done. 
The  camera  was  originally  built  by  R.  Aikens  of  Photometries,  Ltd.,  of  Tucson, 
using  designs  developed  at  Kitt  Peak  National  Observatory. 

Digital  Video  Processor:  This  is  a  partial  duplicate  of  the  SOUP  on-board 
image  processor.  It  receives  digital  frames  from  the  CID  camera,  performs 
arithmetic  for  gain,  offset  and  ratio  corrections  on  three  full-frame 
memories,  and  displays  two  images  simultaneously  on  separate  video  channels. 

A  low-speed  parallel  interface  to  the  microcomputer  is  used  for  commands  and 
(limited)  data  transfer:  passing  a  full  frame  to  the  microcomputer  and 
storing  it  on  disk  takes  about  30  seconds.  Therefore,  videotape  or  videodisk 
are  the  preferred  storage  media. 

Microcomputer  Controller:  The  computer  is  a  DEC  LSI  11-23,  mounted  in  a  VT- 
100  terminal.  It  runs  the  RSX-11M  operating  system  with  a  Winchester  disk  (8 
megabytes)  for  storage.  The  computer  controls  all  hardware  through  DR11C 
parallel  interfaces  to  the  tunable  filter  electronics  box.  Early  versions  of 
the  software  were  written  in  assembler  language.  However,  in  1984,  a  modular 
Fortran  version  was  written  which  controls  the  tunable  filter,  prefilters, 
film  and/or  CID  cameras,  video  processor  and  video  disk  recorder.  New 
scientific  observing  sequences  can  be  entered  at  a  fairly  high  level,  a  vast 
improvement  over  the  previous  assembler  code. 

OBSERVATIONS 

The  first  observing  run  took  place  in  late  October  and  early  November  of 
1982.  Since  a  film  camera  was  used  in  place  of  the  CID  camera,  it  was 
primarily  supported  by  other  funds.  The  LSI-11  computer  orchestrated  the 
tunable  filter,  prefilters  and  35-mm  Acme  camera  in  a  complex  observing 
program  covering  four  solar  lines:  Fe  I  5250  and  5576,  Mg  I  5173  and 
H-alpha.  Two  days  of  very  good  seeing  were  spent  on  a  rapidly  growing  active 
region  near  disk  center.  Approximately  150  of  these  frames  have  been 
digitized  using  the  SPO  fast  microphotometer.  Some  examples  are  included 
in  Appendix  A. 

In  January,  1983,  a  second  observing  run  was  made  to  test  the  Lockheed  active 
mirror  system  of  Dr.  R.  Smithson.  No  spectral  imaging  was  attempted.  We 
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concluded  that  the  active  mirror  was  too  new  and  complex  to  try  a  joint 
observing  run  with  both  instrument  systems  in  1983. 

The  entire  spectral  imaging  system  was  operated  at  the  Tower  Telescope  in 
November,  1983.  An  integrated  control  program  made  a  continuous  series  of 
magnetograms  and  intensity  frames  in  the  blue  wing  of  Fe  I  6302.  These  were 
displayed  on  video  monitors  by  the  image  processor  and  recorded  on  video¬ 
cassettes.  Simultaneously,  Dr.  Luc  Dame  of  LSPS,  France,  recorded  K  line 
images  with  a  Lockheed  filter  and  the  SPO  CCD  camera .  Synchronized  movies 
of  the  K  line  and  magnetograms  with  a  cycle  time  of  12  seconds  were  obtained 
for  about  20  minutes  in  good  seeing.  Later  in  November,  about  30  of  the  best 
magnetograms  were  digitized  from  the  videotape.  A  representative  sample  of 
the  digital  CID  frames  was  also  saved  on  the  Winchester  disk. 

The  fourth  observing  run,  in  December,  1983,  was  a  repeat  of  the  November 
experiment.  Two  hours  of  data  in  fairly  good  seeing  was  obtained.  A  short 
gap  occurred  every  20  minutes  or  so  for  a  tunable  filter  calibration.  The 
video  digitizer  at  Big  Bear  Solar  Observatory  has  been  used  to  obtain  240 
magnetograms  in  digital  FITS  files  on  tape.  These  data  have  not  yet  been 
analyzed;  they  are  being  sent  to  Dr.  Dame  for  power  spectral  analysis  and 
comparison  with  the  K  line  time  series. 

In  September  and  October,  1984,  a  fifth  observing  run  took  place,  supported 
by  other  funds.  This  represented  the  first  serious  attempt  to  operate  the 
tunable  filter  system  in  the  compensated  image  formed  by  the  Lockheed  active 
mirror.  Unfortunately,  this  attempt  was  frustrated  by  bad  weather  in  the 
first  week  of  the  observing  run.  In  the  second  week,  the  complete  spectral 
imaging  system  was  run  in  a  training  exercise  for  the  Spacelab  2  payload 
crew.  The  microcomputer  was  programmed  to  run  15  different  observing 
sequences  simulating  the  SL2  SOUP  instrument.  Only  one  day  of  fair  seeing 
occurred.  The  data  was  videotaped  but  is  not  worth  analyzing  from  a 
scientific  point  of  view. 
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APPENDIX  A 

Examples  of  the  narrowband  solar  Images  obtained  un^er  this  contract  are 
presented  here.  All  of  these  prints  were  copied  from  the  video  of  the 
Grinnell  GMR-274  digital  image  display  at  Lockheed.  A  Tektronix  4634  Image 
Forming  Module  was  used  for  quick,  medium-quality,  gray-scale  prints.  Black 
and  white  glossies,  color  prints  (pseudo-color  renditions)  and  slides  can  also 
be  made  on  request  using  a  Videoprint  5000  camera.  Several  different  routes 
have  led  from  the  telescope  to  the  Grinnell.  Film  frames  (Kodak  2415)  have 
been  developed  at  SPO  and  digitized  with  the  SPO  Fast  Microphotometer  or  the 
Lockheed  PDS  1010G.  A  few  digital  frames  were  stored  on  Winchester  disk  from 
the  video  processor.  Video  magnetograms  recorded  on  Betamax  cassettes  from 
the  video  processor  have  been  digitized  at  Big  Bear  Solar  Observatory.  Mhny 
of  these  images  were  created  in  the  Lockheed  VAX  or  PDP  11  computers  from  raw 
intensity  frames. 
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Magnetogram  of  a  unipolar  plage  in  a  young  active 
region  near  disk  center.  Tick  marks  are  2  arcseconds; 
f ield-of-view  is  approximately  50  x  50  arcseconds. 

The  paper  on  page  B-4  is  based  on  analysis  uf  this 
region,  which  was  observed  in  the  line  Fe  I  5250.2. 
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Continuum  image  (near  5250)  of  the  same  region. 
Three  pores  and  much  abnormal  granulation  are 
present.  All  views  of  this  region  are  cosligned 
to  about  0.25  arcseconds. 
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Velocity  in  the  wings  of  5250,  at  about  +/-  50  mA 
for  the  same  region.  Upflow  is  light,  downfl aw 
Is  dark.  The  zero  point  is  taken  to  be  trc  average 
over  the  whole  frame. 
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Overlay  of  the  locus  of  strong  magnetic  fields  on 
the  Dopplergram  (dark  =  downflow).  Magnetic  flux 
tubes  tend  to  occur  in  downflow  lanes. 
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Overlay  of  the  locus  of  strong  magnetic  fields  on  the 
continuum  image.  The  correlation  with  the  dark  lanes 
is  present  in  some  areas  of  the  picture  but  is  mt 
nearly  as  clear  as  on  the  velocity  image. 
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Mag  netogram  in  Fe  I  6302  made  with  the  CID  camera 
and  video  processor  on  the  November,  1983,  observing 
run.  This  frame  was  digitized  from  videotape. 

It  shows  the  best  seeing  obtained  in  any  run  with 
the  CID  camera.  Tick  marks  are  5  arcseconds. 
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Magnetogram  made  18  minutes  before  the  previous  one. 
Very  few  changes  in  network  fragments  are  seen  at 
this  resolution;  some  of  the  features  inside  the 
large  cell  may  have  changed.  According  to  Dame,  these 
correspond  to  "K-grains"  in  the  Ca  II  K  line. 
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CID  image  of  the  same  region  at  6302.5  line  center. 
This  frame  was  stored  digitally  on  the  Winchester 
disk.  No  gain  correction  has  been  done.  Most  of  this 
fixed  pattern  noise  is  absent  in  the  magnetograms,  but 
it  was  difficult  to  remove  from  the  intensity  frames. 
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Mag  netogram  of  the  same  region,  in  somewhat  worse  seeing 
than  the  previous  ones.  This  frame  was  stored  digitally 
on  the  Winchester  disk.  Variations  in  the  checkerboard 
noise  pattern  of  the  device  cause  the  bad  pixels  on  the 
right  side. 
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Doppler  image  of  the  same  region,  made  at  +/-  50  mA 
in  5250.  Note  the  lack  of  flow  in  the  sunspots, 
except  for  some  Evershed  flow  in  the  penumbrae. 

In  this  region,  the  magnetic-downflow  correlation 
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Continuum  image  of  the  same  region.  This  is  first  of 
a  series  of  frames  showing  the  height  coverage  possible 
with  the  tunable  filter. 
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Image  of  the  same  region  in  the  blue  wing  of  Mg  I  5173, 
showing  the  upper  photosphere  near  the  temperature 
mi n imum. 
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Image  of  the  same  region  in  H-alpha  line  center, 
showing  the  chromosphere.  Note  the  very  bright 
plage  in  the  vicinity  of  the  newly  erupting  bipolar 
magnetic  fields. 


LMSC/F059543 


V 


Magnetogram  of  a  large  cell  on  the  edge  of  an  active 
region,  made  in  Mg  I  5173  in  very  good  seeing.  Tick 
marks  are  2  arcseconds.  This  represents  the  best  seeing 
of  the  1982  observing  run;  many  flux  tubes  can  be  found 
at  the  pixel  limit  of  resolution,  0.5  -  0.75  arcseconds. 
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Blue  wing  intensity  in  Mg  I  5173  of  the  same  region. 

A  time  series  of  9  steps,  spanning  25  minutes,  exists 
but  only  the  first  and  last  step  have  been  digitized. 
Seeing  remains  nearly  this  good  throughout. 
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APPENDIX  B  -  ABSTRACTS 


Abstract  of  Paper  Presented  at  American  Astronomical  Society  Meeting,  Boston, 
January  1983. 


Fine  Structure  in  the  Solar  Magnetic  Field, 1 
T.  D.  Tarbell,  A,  M.  Title,  Lockheed .  We  present  obser¬ 
vations  made  with  the  Spacelab  2  engineering  model  tun¬ 
able  filter  at  the  Sac.  Pk.  Tower  Telescope.  The  filter 
bandwidth  is  75  m£  at  Fe  I  5250,  narrow  enough  to  pro¬ 
vide  some  line  profile  information,  Images  are  recorded 
on  film  and  analyzed  digitally  to  make  longitudinal  mag¬ 
netograms,  Dopplergrams  and  photometrically-corrected 
filtergrams.  The  best  frames  have  resolution  better  than 
0.5  arcseconds.  Results  for  several  active  regions  are 
presented,  illustrating  the  following:  offsets  between 
intensity  and  magnetic  maps  of  faculae,  predicted  by 
magnetostatic  flux  tube  models;  variability  in  the  rela¬ 
tion  between  downdraft  speed  and  magnetic  flux;  the  spa¬ 
tial  relationship  between  granulation  and  flux  tubes  of 
different  sizes;  greater  inclination  of  field  lines  in 
dark  penumbral  filaments  than  in  bright  ones;  extension 
of  the  penumbral  fields  far  beyond  the  visible  penumbral 
edge;  a  possible  change  in  the  magnetic  fine  structure 
of  a  complex  active  region  during  a  flare.  This  work  is 
supported  by  Lockheed  IR  funds  and  AFGL  contracts. 
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Abstract  of  Paper  Presented  at  the  Fourth  European  Meeting  on  Solar  Physics, 
"The  Hydromagnetics  of  the  Sun,”  Utrecht,  October,  1984. 


THE  RELATIONSHIPS  BETWEEN  THE  MAGNETIC  FIELD 
AMD  THE  CHROMOSPHERIC  ‘BRIGHT  POINTS" 

L.  DAME1,  K.  TOPKA2,  A.  TITLE2,  T.  TAjABELL2  and  P.  GOUTTEBROZE1 . 

(1)  LPSP ,  BP  10  ,  Verrieres  le  Buisson  91371,  FRANCE 

(2)  LPARL,  3251  HANOVER  St-,  PALO  ALTO  CA  94304,  USA 

Sirm.il taneous  high  resolution  observations  of  the  magnetic  field  (Lockheed 
high  resolution  Tunable  Lyot  Filter''  and  0-6  A  bandpass  Call  Roy 
Filtergrams  ( l.ocV  h^ed  Biregrlngent  Filter)  were  obtained  in  deeember  1983 
at  the  .'ficcum  Tower  Te'.t-srope  of  Sacramento  Peak  Observatory  on  a  2  by  2 
arcmin  area  of  a;:  lve  and  quiet  sun.  The  excellent  quality  and  high 
resolutions  notn  spatial  ub--arcsec)  and  temporal  (12  sec  exposure  rate) 
on  a  20  minute  sequence  a!  w  a  detailed  comparison  of  the  magnetic  field 
occurence  and  tne  ch r nmn s phc r 1 c  structures  both  in  the  network  and  in  the 
supergranul ar  cells  where  the  magnetic  field  presence  is  compared  to  the 
Call  "Bright  Points"  locations.  The  influence  of  the  chromospheric 
oscillatory  field  on  the  different  structures  and  its  implications  on  the 
correlation  magnetic  field  -  Call  intensities  is  also  investigated. 
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Abstract  of  Paper  Presented  at  American  Astronomical  Society  Meeting,  Tucson, 
January,  1985 

A  mnuscript  describing  this  work  is  being  prepared  for  The  Astrophysical 
Journal. 


High  Resolution  Observations  of  Magnetic  Features 
on  the  Sun 


K.P.  Topka  and  T.D.  Tarbell  (Lockheed) 

The  study  of  small,  isolated,  magnetic  fea¬ 
tures  on  the  sun  reveals  that  they  can  apparent¬ 
ly  decay  in  place,  without  the  simultaneous 
removal  of  an  equal  amount  of  observable  opposite 
polarity  flux.  These  results  were  obtained  from  a 
time  sequence  of  16  high  resolution  magnetograms 
(best  have  1  arc  sec  resolution)  covering  1  hour. 
The  observations,  made  at  Sac  Peak  Observatory, 
made  use  of  the  Lockheed  Tunable  Filter  (engineer¬ 
ing  model  for  the  Spacelab  2  SOUP  experiment  and 
an  active  optical  system  (sunspot  tracker),  during 
moments  of  good  to  excellent  seeing.  Differential 
geometric  image  distortion,  due  to  atmospheric 
seeing,  was  removed  using  a  cross-correlation 
technique.  This  allows  for  a  more  accurate  point- 
by-point  comparison  of  the  magnetograms. 

Twenty-eight  small  (2  to  8  arc  sec  diameters) 
magnetic  features  were  selected  for  study  on  the 
oasis  that  they  were  usually  monopolar  and  well 
isolated  from  surrounding  features.  The  total 
magnetic  flux  present  was  carefully  measured  for 
each  feature  on  each  magnetogram.  About  1/3  of  all 
of  these  features  showed  evidence  for  a  decrease 
in  the  total  amount  of  magnetic  flux  present  dur¬ 
ing  the  1  hour  observing  period,  without  any 
measurable  simultaneous  loss  of  an  equal  amount 
of  nearby  opposite  polarity  flux.  In  one  case  the 
decaying  magnetic  feature  was  very  well  isolated 
from  any  other  feature,  surrounded  by  a  region 
with  no  detectable  magnetic  flux.  During  its 
decay,  no  detectable  magnetic  flux  was  observed 
to  spread  out  into  this  surrounding  region. 
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Abstract  of  Paper  Presented  at  American  Astronomical  Society  Meeting,  Tucson, 
January  1985. 

A  nmnuscript  describing  this  work  is  being  prepared  for  The  Astrophysics  1 
Journal. 


Magnetic  Fields/  Downdrafts,  and  Granulation  in 
the  Solar  Photosphere  "" 

A.M.  Title  and  T.D.  Tarbell  (Lockheed  PARL) 

We  discuss  very  high  resolution  observations  of 
quiet  and  active  sun  in  the  photospheric  Fe  I 
line  5250.  These  images  were  obtained  on  film  at 
the  Sacramento  Peak  Tower  with  a  Lockheed  tunable 
birefringent  filter.  When  magnetograms  and 
dopplergrams  are  coaligned  and  compared  at  this 
resolution  (1/2  arc  second),  we  find  that  the 
weaker  magnetic  flux  concentrations  are  well 
correlated  with  downdrafts.  This  occurs  for 
magnetogram  signals  in  the  100  -  500  gauss  range, 
presumably  corresponding  to  area  filling  factors 
of  10  -  50  percent.  These  fields  form  a  small 
scale  network  pattern,  with  "cell"  diameters  of  2 
to  5  arcseconds;  these  are  somewhat  larger  than 
the  mean  intergranular  spacing  of  2.4  arcseconds. 
Magnetic  fields  in  this  network  correlate  better 
with  dark  intergranular  lanes  than  with  bright 
granules  in  the  coaligned  continuum  images.  This 
relationship  is  less  clear  than  the  downflow 
correlation,  because  our  observed  granule-lane 
contrast  is  much  less  in  magnetic  than  in 
non-magnetic  photosphere  (abnormal  granulation) . 
In  stronger  flux  concentrations  (500  -  1000 
gauss),  the  downflow  is  suppressed  compared  with 
the  weaker  magnetic  surround.  These  correspond  to 
pores  and  "holes"  in  the  continuum  image.  By 
holes,  we  mean  transient  spaces  between  the 
bright  granules  that  are  somewhat  darker  than 
average  but  not  obvious  pores.  The  time-scale  for 
rearrangement  of  the  small  scale  magnetic  network 
is  not  currently  known,  but  it  is  definitely 
longer  than  the  5  to  10  minute  time-scale  of 
granulation. 
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APPENDIX  C 


Paper  Presented  at  the  NSO  Conference,  "Small-Scale 
Dynamical  Processes  in  Quiet  Stellar  Atmospheres," 
Sunspot,  NM,  July,  1983. 


HIGH  SPATIAL  RESOLUTION  MAGNETIC  OBSERVATIONS  OF  AN 
ACTIVE  REGION 


Ken  Topka  and  Ted  Tarbell 

Lockheed  Palo  Alto  Research  Laboratory 

Dept  52-13,  Blgd .  202 

3251  Hanover  St 

Palo  Alto,  California  94304 


ABSTRACT 


High  spatial  resolution  magnetograms  of  an  active 
region  reveal  apparent  changes  in  the  magnetic  flux  of  small 
isolated  features  on  time  scales  less  than  30  minutes.  Both 
flux  increases  and  decreases  are  observed,  sometimes  for  a 
single  feature.  These  changes  are  seen  in  small  unipolar 
magnetic  regions  with  no  obvious  changes  observed  in  any 
nearby  opposite  polarity  features. 


1.  INTRODUCTION 

The  birth  of  an  active  region  always  results  from  the 
emergence  of  new  magnetic  flux  from  below  the  photosphere, 
never  from  the  concentration  of  previously  existing 
photospheric  flux.  Many  new  active  regions  are  born  each 
sunspot  cycle,  and  this  implies  that  a  lot  of  new  flux  must 
emerge.  An  equal  amount  of  flux  must  be  removed  from  the 
photosphere  each  cycle  as  well,  in  order  to  prevent  an  ever 
increasing  buildup  of  flux  on  the  sun. 

It  has  been  known  for  some  time  that  active  regions 
decay  on  time  scales  of  a  few  months.  The  flux  appears  to 
spread  out  by  diffusion  (cancelling  where  opposite 
polarities  meet) ,  until  that  which  remains  is 
indistinguishable  from  the  background.  More  recently, 
significant  flux  changes  in  active  regions  on  time  scales  of 


a  few  days  have  been  reported  (Wallenhorst  and  Howard  1981) 
(Wallenhorst  and  Topka  1982).  Wilson  and  Simon  (1983)  have 
reported  flux  changes  in  small  magnetic  regions  on  time 
scales  of  about  20  minutes. 

The  purpose  of  this  paper  is  to  report  on  the  search 
for  flux  changes  in  small  magnetic  regions  on  time  scales  of 
a  few  minutes,  similar  to  those  reported  by  Wilson  and 
Simon.  The  observations  analyzed  were  selected  on  the  basis 
of  good  time  coverage  and  good  seeing,  rather  than  choosing 
a  specific  target  or  event. 


2.  OBSERVATIONS 

The  observations  were  obtained  on  19  April  1980  (1427 
UT  -  1453  UT)  of  AR  2391,  an  older  active  region  located  at 
N  12,  W  44 .  The  equipment  consisted  of  the  Lockheed  tunable 
filter  (100  mA  bandpass  FWHM,  engineering  model  for  Spacelab 
2  Solar  Optical  Universal  Polarimeter;  Title  and  Rosenberg 
1981) ,  mounted  behind  an  image  stabilization  system  that 
guides  on  a  sunspot  (Tarbell  and  Smithson  1981) .  The 
observations  were  carried  out  at  the  Sacramento  Peak 
Observatory  Vacuum  Tower  Telescope. 

All  images  were  recorded  on  photographic  film. 
Fourteen  images  were  recorded  at  each  time  step,  six  in  the 
magnetically  sensitive  Fe  I  6302.50  A  line  (3  RCP  +  LCP 
pairs) •  and  eight  in  H-alpha  (wavelength  scan  from  H-alpha  + 
2  A  to  H-alpha  -  2  A) .  Magnetic  pairs  were  taken  in  the 
blue  wing  of  6302,  60  mA  from  line  center.  The  filter 
wavelength  was  calibrated  periodically  with  a  HeNe  laser  and 
found  to  be  stable  to  within  a  few  mA.  The  seeing  was  quite 
good  and  fairly  uniform,  permitting  one-half  to  one 
arcsecond  resolution  despite  the  relatively  long  exposure 
times  needed  (0.5  seconds). 


3.  RESULTS 

Fifty  four  images  comprising  9  time  steps  (  3  Fe  I  RCP 
+  LCP  nairs  per  time  step)  have  been  digitized  (256  x  256 
pixels)  using  the  fast  microphotometer  at  Sac  Peak.  Film 
densities  are  converted  to  intensities  using  a 
characteristic  curve  derived  from  exposure  tests.  One 
digital  image  at  each  time  step  was  arbitrarily  chosen  as  a 
standard.  Each  of  the  other  five  images  at  the  same  time 
step  were  subdivided  into  16  equal  parts,  and  each  part  was 
matched  as  closely  as  possible  (via  cross-correlation)  to 
the  corresponding  n*rt  of  the  standard  image.  After  this 
was  performed  for.  all  16  parts  of  each  image,  they  were 
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stretched  (bilinear  interpolation  between  the  measured 
offsets  in  X  and  Y)  to  correspond  as  closely  as  possible  to 
the  standard.  This  process  does  not  remove  all  image 
distortion  due  to  seeing,  instead  it  stretches  5  of  the 
images  until  they  closely  match  the  image  distortion  present 
in  the  sixth  image  at  the  same  time  step.  The  match  is 
generally  good  to  better  than  1/4  to  1/2  pixel.  Each  RCP  + 
LCP  pair  has  been  subtracted  to  produce  27  magnetograms. 
Point  by  point  subtraction  of  a  pair  of  images  (1  RCP,  1 
LCP)  leads  to  the  determination  at  each  point  of  the  V 
Stokes  parameter  (  V  =  [RCP  -  LCP]).  For  a  uniform  magnetic 
field  parallel  to  the  line  of  sight  the  magnetic  field 
strength  B  is  approximately  proportional  to  V.  More 
generally,  the  normalized  magnetogram  V/I  (where  I  =  RCP  + 
LCP)  is  approximately  proportional  to  the  magnetic  flux 
averaged  over  the  resolution  element  of  the  image.  The  3 
magnetograms  at  each  time  step  were  then  averaged  to  produce 
9  final  magnetograms,  one  at  each  time  step.  The  field  of 
view  is  80  x  80  arc  sec,  and  the  limiting  spatial  resolution 
is  about  1  arc  sec. 

Twenty  two  small,  unipolar  magnetic  features  were 
selected  for  study.  Most  features  were  isolated,  to  avoid 
apparent  flux  changes  that  are  due  to  undetectable  movement 
of  magnetic  features  into  and  out  of  the  region  of  study. 
Six  were  selected  on  the  basis  of  visual  inspection  of  the 
magnetograms  (they  appeared  to  change) ,  and  the  rest  were 
picked  at  random.  A  rectangular  box  was  drawn  around  each 
feature,  and  the  average  Stokes  V/I  was  determined  in  each 
box  for  all  9  magnetogram  frames.  The  average  Stokes  V/I 
was  multiplied  by  the  area  of  each  box  to  yield  the  net 
magnetic  flux  present  (in  arbitrary  units,  since  all  we  are 
interested  in  here  are  flux  changes.  Later  we  discuss 
estimates  of  the  actual  amount  of  magnetic  flux  involved) . 
The  boxes  ranged  in  area  from  4  to  40  square  arc  sec, 
depending  upon  the  size  of  the  magnetic  feature  enclosed. 
The  boxes  were  made  large  enough  so  that  slight  shifts  of  a 
magnetic  feature  (from  one  magnetic  frame  to  the  next)  due 
to  image  distortion  induced  by  atmosphere  seeing  does  not 
seriously  affect  the  results. 

There  were  significant  fluctuations  in  the  zero  points 
of  the  magnetograms  from  one  frame  to  the  next,  due  in  part 
to  slight  changes  in  density  with  position  on  each  original 
photographic  frame.  A  local  background  (a  field  free  area 
that  is  as  close  as  practical  to  each  feature)  was  chosen. 
Subtraction  of  the  local  background  significantly  reduced 
frame  to  frame  fluctuations,  and  the  typical  small  feature 
showed  a  maximum  variation  of  about  10%  in  the  measured  flux 
values,  with  no  trends  discernable.  Four  small  magnetic 
regions  with  no  flux  changes  greater  than  the  fluctuations 
described  above  are  shown  in  figure  1. 


Other  features,  however,  showed  systematic  changes  in 
the  measured  magnetic  flux  that  were  much  greater  then  the 
10%  maximum  found  for  most  small  magnetic  regions.  One 
feature,  for  example,  lost  almost  90%  of  its  flux  in  about 
10  minutes.  Figure  2  shows  5  examples  of  features  with 
apparent  changes  in  flux.  Included  in  figure  2  are  flux 
decreases,  increases,  and  both  increases  and  decreases  for  a 
single  feature  (flux  decay  followed  by  recovery) .  In  one 
case  the  possible  movement  of  flux  from  outside  to  inside  of 
the  box  containing  a  feature  may  explain  some  of  the 
observed  flux  change.  Feature  5  appeared  to  decay  while 
new  flux  emerged  nearby.  Some  of  this  new  flux  appears  to 
have  migrated  into  the  area  near  feature  5,  thus  causing  an 
apparent  recovery  of  flux  (see  figure  2) .  The  box 
containing  feature  5  was  only  about  4  square  arc  seconds  in 
area,  and  it  was  not  possible  to  shift  it  around  to 
completely  avoid  the  new  emergence  without  excluding  part  of 
feature  5.  No  evidence  for  flux  migration  was  seen  for  the 
other  features.  The  5  features  of  figure  2  had  the  most 
pronounced  trends  of  the  22  studied,  and  all  were  selected 
because  they  appeared  to  be  changing  visually.  Most  of  the 
rest  of  the  features  studied  had  no  discernable  trends  in 
their  measured  magnetic  fluxes,  similar  to  the  4  features 
shown  in  figure  1.  Less  than  5%  of  the  total  magnetic  flux 
present  within  our  magnetogram  field  of  view  is  contained  in 
these  22  features. 

Table  1  lists  estimates  of  the  amount  of  magnetic  flux 
initially  present  in  each  of  the  features  shown  in  figure  2, 
as  well  as  the  change  in  flux,  the  time  for  this  change  to 
occur,  and  the  rate  of  change  of  flux.  The  magnetogram 
signals  (V/I)  have  been  converted  to  magnetic  flux  using  the 
slope  of  the  6302  line  profile  in  similar  magnetic  elements, 
measured  from  wavelength  scans  with  the  same  filter.  While 
this  calibration  could  be  wrong  by  as  much  as  50%,  we  are 
concerned  here  with  the  temporal  changes,  not  with  the 
absolute  values.  From  the  line  profile  scans,  we  can  infer 
that  changes  in  Doppler  shifts  of  several  km/sec  would  be 
needed  to  affect  the  flux  measurements  at  the  levels 
observed. 
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Figure  1.  The  magnetic  flux  (or  Stokes  parameter  V  for 
the  upper  left  plot)  is  shown  as  a  function  of  time  for 
4  different  isolated  features.  The  straight  lines 
through  the  data  points  represent  the  average  value. 
No  strong  flux  increases  or  decreases  are  apparent  in 
these  data  . 
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Figure  2.  The  magnetic  flux  for  5  different  isolated 
features  is  shown  plotted  as  a  function  of  time.  Both 
flux  increases  and  decreases  as  apparent,  as  indicated 
by  the  trend  lines. 


TABLE  1.  MEASURED  FLUX  CHANGES  IN  SMALL  MAGNETIC  REGIONS 


FEATURE  INITIAL  FLUX  CHANGE  IN  FLUX  TIhE  RATE  OF  CHANGE 
(x  1.0E+18  Mx)  (x  1.0E+18  Mx)  (min)  (1.0E+15  Mx/s) 


10.2 


10.0 


-5.4 


-3  .7 


+  2  .7 


+7.7 


+  9.2 
-3.3 


18.8 


26.0 


26.0 


-8.6 


+2.4 


+4.9 


+5.9 


-5.7 


Note:  The  initial  flux  given  in  column  2  is  absolute  value. 

The  -  sign  in  column  3  denotes  flux  decreases,  +  sign 
flux  increases. 


4.  DISCUSSION 

Because  the  divergence  of  the  magnetic  field  vector  B 
is  zero,  an  eruption  of  positive  magnetic  polarity  flux  must 
be  accompanied  by  a  simultaneous  and  equal  eruption  of 
negative  polarity  flux.  The  typical  picture  for  this 
process  is  the  emergence  through  the  photosphere  of  a  loop 
consisting  of  one  or  more  'flux  tubes',  due  to  the  magnetic 
bouyancy  of  the  loop  and/or  convective  upflows.  As  the  loop 
emerges  two  footpoints  of  opposite  polarity  form  that  move 
apart  (Zwaan  1978).  This  process  has  been  observed  in 
similar  filter  data  of  other  active  regions. 


The  long  term  decay  of  an  active 
to  diffusion.  Since  magnetic  field 
the  photospheric  and  subphotospher ic 
carried  along  with  any  flows  in  this 
field  lines  are  assumed  to  random- 
points  of  emergence,  and  magnetic 
spread  out  due  to  Ohmic  dissipat 
polarities  eventually  meet  they  ' 
results  in  the  net  removal  of  flux  fr 


region  is  probably  due 
lines  are  'frozen'  into 
material,  they  will  be 
material.  The  magnetic 
walk  away  from  their 
features  may  further 
ion.  When  opposite 
cancel',  a  process  that 
om  the  photosphere. 


However,  there  now  exist  observations  that  suggest  that 
processess  other  than  those  described  above  are  present  on 


the  sun.  The  process  of  emergence  has  been  observed  where 
only  one  polarity  flux  is  observable,  thus  contradicting  the 
view  the  two  footpoints  of  equal  and  opposite  polarity  that 
move  apart  must  always  be  seen.  In  some  cases  no  observable 
increase  in  the  flux  of  the  opposite  polarity  is  seen 
anywhere  near  the  site  of  the  emergence  (and  in  fact,  in 
some  cases  no  opposite  polarity  flux  is  detectable  at  the 
site  of  the  emergence  at  all) .  The  observations  of  Wilson 
and  Simon  (1983)  include  this  phenomenon,  as  do  the 
observations  reported  in  this  paper.  We  do  not  suggest  that 
magnetic  monopoles  exist  on  the  surface  of  the  sun.  The 
missing  polarity  flux  is  there,  somewhere,  and  may  be 
undetectable  in  our  magnetograms  for  some  reason.  In  any 
case,  the  phenomenon  is  not  yet  understood. 

The  decay  of  an  active  region  via  diffusion  requires 
long  time  scales  (months)  for  the  removal  of  significant 
amounts  of  photospheric  flux.  Recent  observations  show 
possible  flux  removal  from  active  regions  on  much  shorter 
time  scales  than  is  expected  from  diffusion,  again 
suggesting  other  mechanisms  exist.  Howard  and  La  Bonte 
(1981)  made  the  rather  startling  claim  that  fully  10%  of  the 
magnetic  flux  present  on  the  sun  today  was  not  there 
yesterday,  implying  that  emergence  and  disappearance  of  flux 
to  and  from  the  photosphere  is  nearly  continuous,  and  that 
large  amounts  can  be  added  or  subtracted  on  relatively  short 
time  scales.  Wallenhorst  and  Howard  (1982)  and  Wallenhorst 
and  Topka  (1982)  observed  that  a  substantial  fraction  of  the 
entire  flux  of  an  active  region  can  apparently  be  removed 
from  the  photosphere  in  just  a  few  days.  The  observed  decay 
rates  were  up  to  10%  of  the  total  flux  of  the  active  region 
per  day.  Furthermore,  Wallenhorst  and  Topka  (1982)  obtained 
higher  resolution  magnetograms  (best  were  2  arc  sec,  typical 
was  2-3  arc  sec)  .  These  magnetograms  seemed  to  show  a 
general  weakening  of  the  active  region,  even  in  large 
unipolar  areas  far  removed  from  the  neutral  line.  Magnetic 
flux  was  just  fading  away  in  place,  rather  than  coming 
together  with  equal  and  opposite  polarity  flux  at  the 
neutral  line  and  disappearing  there. 

Wilson  and  Simon  (1983)  reported  flux  decreases  on  very 
short  time  scales  (less  than  1  hour),  without  any  observable 
changes  in  nearby  opposite  polarity  flux.  In  these 
observations  the  entire  line  profile  (6302  A)  was  used, 
instead  of  just  a  narrow  band  pass  in  the  wings  of  the  line 
(as  in  a  conventional  Babcock  magnetograph) .  It  is  hard  to 
see  how  instrumental  effects  or  line  profile  effects  could 
have  caused  the  flux  changes.  Although  the  observations 
reported  here  are  of  the  magnetograph  type  that  uses  only 
the  blue  wing  of  the  line,  they  show  essentially  the  same 
sorts  of  flux  increase  and  decay.  The  phenomena  of 
emergence  and  of  'disappearance  in  place'  of  one  polarity  of 
flux,  without  any  observable  opposite  polarity  flux  emerging 


or  disappeari  j  nearby,  has  potentially  interesting 
theoretical  consequences  if  confirmed. 
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ABSTRACT 

We  present  a  critical  analysis  of  observations  of  the  decay  of 
large  active  regions/  and  we  propose  that  magnetic  flux  disappears 
more  rapidly  than  can  be  explained  with  Leighton's  diffusion  model. 
The  proposed  component  of  the  flux  disappearance  seems  to  take  place 
locally,  without  any  apparent  motion  of  magnetic  elements,  and  does 
not  seem  to  involve  the  submergence  of  magnetic  flux. 

To  explain  these  observations  we  propose  a  model,  in  which  part 

of  the  observed  large-scale  field  of  active  regions  is  annihilated  by 

small-scale  fields  that  emerge  in  between  patches  of  opposite 

polarity.  The  small-scale  field  is  assumed  to  consist  of  a  train  of 

17 

bipoles,  each  with  typical  size  of  1000  km  and  flux  of  10  Mx,  and 
with  an  orientation  opposite  to  that  of  the  active  regions  as  a  whole. 
The  interaction  of  these  bipoles  with  the  overlying  coronal  field 
leads  to  the  decay  of  the  active  region,  without  any  apparent 
large-scale  motion  of  magnetic  features  in  the  photosphere. 

To  show  that  such  a  scenario  is  not  physically  unreasonable,  we 
describe  a  possible  mechanism  by  which  the  small-scale  field  might  be 
formed.  According  to  this  model,  part  of  the  active-region  field 
reconnects  deep  below  the  solar  surface,  leading  to  the  formation  of 
U-shaped  magnetic  loops  in  the  convective  zone.  During  the  decay  of 
the  active  region,  these  U-loops  float  to  the  surface,  where  they 
emerge  as  trains  of  small  bipoles.  Due  to  the  large  amount  of  mass 
captured  in  the  U-loop,  the  predicted  f ieldstrength  of  the  bi poles  is 
only  a  few  Gauss. 
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Subject  Headings:  hydromagnetics  —  sun:  activity 
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1.  INTRODUCTION 


Magnetic  flux  first  appears  at  the  solar  surface  in  the  form  of 

easily  identified,  bipolar  active  regions.  The  emergence  of  flux  in 

both  large  and  small  active  regions  takes  place  in  a  similar  way,  and 

is  fairly  well  understood  observationally :  a  loop  of  magnetic  field 

breaks  through  the  solar  surface  from  below,  and  the  opposite  polarity 

footpoints  of  the  loop,  as  seen  in  photospheric  magnetograms,  separate 

as  the  loop  rises  up  through  the  atmosphere  (Bumba  and  Howard  1965; 

Zwaan  1978;  Sheeley  1981) .  Presumably,  the  flux  of  active  regions 

originates  from  a  horizontal  field  located  at  some  level  in  the 

1 1 

convective  zone.  Large  active  regions,  with  flux  greater  than  10 
Mx,  probably  originate  from  the  deep  layers  of  the  convective  zone 
(e.g.  Parker  1975.  1984a;  Golub  et  al.  1981). 

In  comparison  with  flux  emergence,  the  disappearance  of  magnetic 
flux  from  the  surface  is  much  less  well  understood  (Zwaan  1978) .  As 
the  field  of  an  active  region  is  dispersed  across  the  solar  surface, 
the  measured  flux  is  seen  to  decrease,  and  ultimately  the  region 
disappears  into  the  background  field.  It  has  been  presumed  that  most 
of  the  flux  disappearance  takes  place  in  regions  where  opposite 
polarity  fields  meet  at  neutral  lines,  and  in  areas  of  mixed  polarity. 
The  fact  that  flux-disappearance  is  not  a  dramatic  process  indicates 
that  it  probably  occurs  on  small  length-scales,  and  involves  only 
small  amounts  of  flux  at  a  time.  Leighton  (1964)  proposed  that  the 
intermixing  of  opposite  polarity  fields,  due  to  random  motions  induced 
by  supergranulation,  plays  an  important  role  during  the  late  phases  of 
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evolution  of  active  regions.  In  Leighton's  model  the  action  of  the 
supergranular  motions  on  the  field  is  described  as  diffusion,  which 
implies  a  certain  relation  between  dispersal  and  decay  of  the  field. 

However,  it  has  been  suggested,  e.g.  by  Wallenhorst  and  Topka 
(1982) ,  that  another  process  of  flux-disappearance  takes  place  as 
well.  The  observations  of  Wallenhorst  and  Topka  suggested  that  part 
of  the  active- region  flux  disappears  locally,  without  any  apparent 
motion  of  field  towards  the  neutral  line;  the  two  patches  of  opposite 
polarity  just  fade  away.  This  seems  in  contradiction  with  our 
theoretical  understanding:  due  to  the  high  conductivity  of  the  solar 
plasma,  the  field  is  expected  to  be  "frozen"  into  the  gas,  so  that  the 
evolution  of  the  magnetic  field  is  due  primarily  to  the  motions  of  the 
plasma;  magnetic  diffusion  is  thought  to  play  an  insignificant  role. 
Then  it  seems  that  the  only  way  in  which  flux  can  be  annihilated  is 
when  elements  of  opposite  polarity  move  towards  each  other  until  they 
merge  and  disappear.  Since  the  opposite  polarities  in  a  large  active 
region  typically  are  50,000  km  away  from  the  neutral  line,  we  would 
expect  to  see  motions  of  magnetic  features  over  at  least  this  distance 
during  the  decay  of  the  region.  The  observed  decay  of  the  flux 
without  any  significant  spreading  or  dispersal  of  the  field  thus  poses 
a  difficult  theoretical  problem,  which  we  will  address  in  the  present 
paper. 

The  paper  is  set  up  as  follows.  In  Sect.  2  we  discuss  various 
observations  of  active- region  decay,  with  the  purpose  to  extract 
further  evidence  for  the  proposition  that  part  of  the  flux  disappears 
locally.  In  an  attempt  to  explain  these  observations  we  propose,  in 
Sect.  3,  a  qualitative  model  of  the  decay  of  a  large  active  region,  in 
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particular  regarding  the  dynamics  of  the  suh-surface  magnetic  field. 
We  consider  an  isolated  active  region,  i.e. ,  we  neglect  possible 
interactions  with  other  active  regions,  ephemeral  regions,  or  other 

background  fields. 


2.  OBSERVATIONS  OF  FLUX  DISAPPEARANCE 


In  the  following  we  discuss  various  observations  related  to  the 
disappearance  of  flux,  in  order  to  identify  the  mechanisms  by  which 
this  process  takes  place. 


a)  Large-scale  characteristics 


In  the  course  of  the  solar  cycle  a  large  number  of  active  regions 

appear,  yet  there  is  no  large  build-up  of  field  at  the  surface.  This 

implies  that  most  magnetic  flux  disappears  on  a  short  timescale.  The 

average  lifetime  of  magnetic  flux  at  the  surface  may  be  estimated  in  a 

number  of  ways.  One  method  is  to  divide  the  total  amount  of  flux 

present  on  the  surface  at  any  one  time  by  the  rate  of  emergence  of 

large  active  regions.  Observations  by  J.W.  Harvey  (1983,  private 

2S 

communication)  indicate  that  on  average  about  4  10  Mx  of  each 

polarity  was  present  on  the  solar  surface  during  the  maximum  of  cycle 

21.  (This  flux-value  has  not  been  corrected  for  the  effects  of 

line-profile  changes  in  sunspots,  nor  for  the  fact  that,  at  high 

latitudes,  only  the  line-of-sight  component  of  the  field  was  observed; 

however,  these  corrections  are  probably  small) .  The  rate  of  flux 

emergence  during  the  Skylab  period  in  1973  was  discussed  by  Golub  et 

al.  (1976,  1979).  They  estimated  that  active  regions  with  lifetimes 

21 

of  a  few  days  or  more  contributed  about  8  10  Mx/day.  For  the  year 
1956,  Sheeley  (1981)  estimated  an  emergence  rate  of  3  10  Mx/day  for 
regions  with  lifetimes  of  at  least  10  days.  These  observations 
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suggest  that  an  emergence  rate  between  about  3  10  and  10  Mx/day  is 
typical  for  solar  maximum.  Combining  this  with  Harvey's  value  for  the 
surface  flux,  we  find  that  the  average  lifetime  of  the  flux  in  large 
active  regions  is  between  1  and  4  months. 

A  somewhat  different  method  to  estimate  the  lifetime  is  that  of 

O 

Howard  and  LaBonte  (1981),  who  measured  flux  changes  in  regions  13.2 

o  V 

in  longitude  and  3.4  /cos  A  in  latitude  over  a  period  of  3  days  near 
central  meridian  passage.  They  found  that  flux  appears  and  disappears 
at  a  rate  that  is  sufficient  to  replace  all  existing  flux  in  about  10 
days,  from  which  we  infer  that  the  flux  typically  changed  by  30%  over 
the  3-day  observing  period.  Positive  and  negative  flux  changes  occur 
in  equal  numbers,  and  with  the  same  amplitude.  Although  this  10-day 
timescale  is  significantly  shorter  than  that  estimated  from  the 
emergence- rate  of  active  regions,  it  should  be  kept  in  mind  that  the 
method  of  Howard  and  LaBonte  (1981)  provides  only  a  lower  limit  to  the 
average  lifetime.  Since  it  is  not  known  how  long  a  period  of  flux 
increase  or  decrease  lasts,  it  is  possible  that  a  10%  flux  decrease 
per  day  persists  only  for  about  3  days,  and  is  followed  by  a  flux 
increase.  Then  the  net  effect  is  a  much  slower  change,  with  a  time 
constant  longer  than  10  days. 

Rapid  disappearance  of  flux  is  also  suggested  by  the  observations 
of  Gaizauskas  et  al.  (1983) ,  who  looked  for  large-scale  patterns  in 
the  emergence  of  solar  active  regions.  Using  synoptic  maps  of  the 
magnetic  field  over  27  solar  rotations  obtained  from  Kitt  Peak 
magnetograms,  they  identified  complexes  of  activity.  These  are  series 
of  apparently  related  active  regions  that  emerge  near  the  same 
longitude  and  latitude  over  a  period  of  3  to  6  rotations.  An 
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important  aspect  of  these  activity  complexes  is  that  they  are 

maintained  by  fresh  injections  of  magnetic  flux  in  the  form  of 

distinct  bipolar  regions,  while  the  fiux  of  earlier  regions 

disappears.  In  the  case  of  one  particularly  large  complex  it  was 

11 

shown  that  the  magnetic  flux  remained  nearly  constant  at  7  10  Mx 

during  6  rotations,  even  though  large  new  regions  kept  appearing  at  an 

average  rate  of  one  every  7  days.  As  Gaizauskas  et  al.  (1983)  point 

out,  this  indicates  that  flux  disappearance  in  this  complex  takes 

place  in  situ,  i.e. ,  the  flux  is  not  dispersed  across  the  solar 

surface.  When  a  complex  dies  the  flux  disappears  on  a  timescale  of 

only  one  or  two  months.  This  also  seems  consistent  with  the 

observations  of  J.W.  Harvey  (1983,  private  communication),  who  found 

22 

that  the  total  flux  of  the  sun  fluctuates  by  about  20%,  or  7  10 
Mx,  over  a  period  of  2  solar  rotations;  these  variations  may  well  be 
due  to  the  birth  and  decay  of  activity  complexes. 


b)  Diffusion  models 


The  spreading  of  active  region  magnetic  fields  across  the  solar 
surface  has  been  studied  by  a  number  of  authors.  Leighton  (1964) 
showed  that  the  dispersal  of  unipolar  and  bipolar  magnetic  regions  can 
be  described  fairly  well  as  a  diffusion  process,  i.e.,  the  flux  h(jl,  t) 
per  unit  area  satisfies  the  diffusion  equation: 
z 

~*7  =  d  V n_  ,  (l) 

where  D  is  the  diffusion  constant,  and  V  is  the  Laplacian  operator  on 
the  surface  of  a  sphere.  Leighton  ascribed  the  diffusive  action  to 
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supergranulation,  which  moves  the  magnetic  elements  about  in  a 
random-walk  fashion.  The  interdiffusion  of  opposite  polarity  field 
brings  elements  of  positive  and  negative  flux  close  together  and  leads 
to  the  disappearance  of  these  elements.  The  rate  at  which  opposite 
polarity  fields  interdiffuse  depends  on  the  magnitude  of  the  diffusion 
constant.  The  value  of  the  diffusion  coefficient  D  used  by  Leighton 
was  between  780  and  1550  km  /s. 

A  more  refined  estimate  of  the  diffusion  constant  was  given  by 
Mosher  (1977)  ,  who  used  several  independent  methods  to  determine  D. 
From  the  rate  of  separation  of  opposite  polarities  in  long-lived 

i 

active  regions  Mosher  estimated  that  D  200  km  /s;  this  estimate 
takes  into  account  the  effect  of  differential  rotation  on  the 
evolution  of  the  flux  distribution  (Leighton,  1964) .  Observations  of 
the  rate  of  area  increase  of  magnetic  regions  are  consistent  with  this 
estimate.  The  observed  lifetimes  of  active  regions  yield  D  =  260  - 
770  kmVs,  depending  on  the  assumed  value  of  the  threshold 
f ieldstrength  below  which  weak  fields  become  undetectable.  Mosher 

(1977)  also  measured  the  diffusion  coefficient  from  the  observed 
motion  of  magnetic  features  over  short  periods  of  time  (hours  to 
days) .  Using  different  data  sets  and  methods  of  analysis,  values  of  D 
between  140  and  300  km  /s  were  found.  Mosher  concluded  that,  within  a 
factor  2,  the  diffusion  constant  D  =  200  km  /s,  and  that  the  spreading 
of  magnetic  field  over  periods  of  months  can  indeed  be  represented  as 
a  diffusion  process. 

However,  Mosher  (1977)  points  out  that  the  decay  of  the  active 

region  flux  seems  to  be  more  rapid  than  can  be  explained  with 

x 

intermixing  due  to  diffusion  with  D  =  200  km  /s.  To  illustrate  his 
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concern,  Mosher  considered  a  simple  model  of  an  active  region,  in 
which  all  positive  flux  is  initially  located  at  a  point  x  =  +a,  and 
all  negative  flux  at  x  =  -a,  on  the  x-axis;  NQ  is  the  initial  source 
strength  in  Maxwells.  The  dispersal  and  intermixing  of  flux  as 
function  of  time  t  is  described  by  Eq.(l),  which  yields  the  solution: 

/v  X  -*.y  t  yo 

n<x,y,t)  "nttr  expI'(~^£r_)1  sinh(iDt)-  (2) 

Here  x  and  y  are  horizontal  coordinates,  and  effects  due  to  the 
curvature  of  the  solar  surface  have  been  neglected.  The  net  surviving 
flux  (of  one  sign)  at  time  t  is  obtained  by  integrating  n(x,y,t)  over 
the  half-plane  x  >  0: 


N+(t>  *N»erf(W’' 


(3) 


where  erf  is  the  error  function  (see  Abramowitz  and  Stegun  1965) . 

According  to  this  model  approximately  half  the  initial  flux  has 

x 

disappeared  at  time  t  =  t  -  a  /D.  However,  the  remaining  flux 

*/i 

disappears  rather  slowly,  since  N+ w  0.56  N0(t0/t)  for  t»to;  for 

the  flux  to  drop  below  20%  of  its  initial  value  we  must  wait  until  t  = 

* 

7.8  t0.  For  an  initial  separation  of  2  a  =  6  10  km  and  a  diffusion 

z 

constant  D  =  200  km  /s,  the  timescale  tQ  =  52  days,  so  that  20%  of  the 
flux  still  remains  after  13  months.  But  in  reality  very  little  flux 
remains,  even  of  large  active  regions,  after  only  6  months. 


This  enhanced  decay  rate  explains  why  the  diffusion  constant, 
deduced  by  Mosher  (1977)  from  the  observed  lifetimes  of  active 

<Z 

regions,  was  somewhat  larger  than  200  km  /s:  to  reproduce  the  actual 
flux-decay  rate  with  a  diffusion  model,  a  larger  diffusion  constant  is 

S3 
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necessary.  Mosher  suggested,  however,  that  no  great  significance 
should  be  attached  to  these  facts,  since  his  flux  estimates  were 
obtained  by  visual  inspection  of  magnetograros  with  a  lowest  contour 
level  of  5  Gauss,  and  consequently  his  flux  measurements  are  rather 
crude.  He  also  proposed  the  existence  of  a  poleward  meridional  flow 
that  transports  the  weak  remnant  fields  out  to  higher  latitudes.  This 
meridional  flow  sweeps  clean  the  low-latitude  active  belts,  thus 
preventing  the  build-up  of  magnetic  fields  there. 

Although  the  presence  of  a  meridional  flow  has  been  reported  by 
other  observers  (Beckers  1979;  Howard  1979;  Duvall  1979;  LaBonte 
and  Howard  1982;  Topka  et  al.  1982),  we  suggest  that  meridional 
motions  alone  are  not  sufficient  to  explain  the  rapid  decay  of  active 
regions.  Also,  we  believe  that  the  difference  between  the  observed 
flux  and  the  flux  predicted  by  the  diffusion  model  (with  D  =  200 
km  /s)  is  too  large  to  be  due  to  observational  uncertainties.  In  our 
opinion  the  lack  of  a  long-term  residue  of  active  regions  points  to  a 
more  fundamental  problem  with  the  diffusion  model. 

Recently,  Sheeley  et  al .  (1983),  Boris  et  al .  (1983),  and  Sheeley 
et  al.(1984,  in  preparation)  have  studied  the  decay  of  a  number  of 
isolated  active  regions,  using  Kitt  Peak  magnetograms.  They 
determined  the  value  of  the  diffusion  constant  D  by  comparing  two 
magnetograms,  taken  one  or  more  rotations  apart,  with  the  results  of  a 
diffusion  model.  This  model  includes  the  effects  of  differential 
rotation  and  meridional  flow  on  the  flux  distribution.  In  contrast 
with  Mosher's  study,  the  method  is  sensitive  both  to  the  shape  of  the 

magnetic  contours  and  the  amount  of  flux  in  the  region.  In  15  out  of 

2. 

18  cases  the  authors  find  values  of  D  between  280  and  600  km  /s,  with 
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an  average  of  430  km  /s,  somewhat  higher  than  the  value  found  by 
Mosher  (1977)  .  The  remaining  three  cases  have  D  >  970  kmVs,  the 
cause  of  which  is  unknown. 

We  suggest  that  the  difference  between  the  430  kmVs  of  Sheeley 

et  al.  (1984)  ,  and  the  200  km^/s  of  Mosher  (1977),  is  real  and  is  due 

to  the  fact  that  Mosher  did  not  try  to  fit  the  observed  flux. 

Although  the  dispersal  of  magnetic  fields  on  the  surface  is  well 

represented  by  diffusion  with  D  =  200  km2/s,  the  rate  of  flux 

disappearance  is  too  large  to  be  explained  in  this  way.  We  suggest  t 

that  a  significant  fraction  of  the  flux  disappears  locally,  i.e., 

without  the  intermixing  of  opposite  polarity  fields.  Fitting  the 

magnetograph  data  with  a  diffusion  model,  Sheeley  et  al.  have 

z. 

determined  an  effective  "diffusion"  constant  (De^  =  430  km  /s)  that 
includes  the  effects  of  enhanced  flux  decay. 

Marsh  (1978)  has  suggested  that  ephemeral  active  regions  (cf. 
Martin  and  Harvey  1979)  may  cause  an  enhancement  of  the  diffusion 
constant.  Interactions  between  ephemeral  regions  and  network  elements 
in  the  quiet  sun  have  recently  been  observed  by  Martin  (1983) . 
Reconnection  of  one  of  the  poles  of  an  ephemeral  region  with  a 
neighboring  network  element  is  equivalent  to  the  displacement  of  this 
element  over  a  distance  on  the  order  of  the  pole  separation  of  the 
region.  Marsh  (1978)  estimated  that  the  effective  diffusion  constant 

X 

due  to  this  process  is  about  800  km  /s.  We  suggest,  however,  that 
ephemeral  regions  are  important  only  for  the  diffusion  of  weak  fields, 
i.e.,  during  the  late  decay-phase  of  an  active  region,  and  in  the 
"quiet"  sun  outside  active  regions.  The  reason  is  that,  in  strong 
magnetic  fields,  ephemeral  regions  do  not  have  sufficient  flux  to 
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affect  the  network  elements.  Marsh  (1978)  made  two  assumptions,  which 
restrict  his  analysis  to  a  particular  f ieldstrength :  (1)  it  was 
assumed  that  ephemeral  regions  fully  reconnect  with  network  elements, 
so  that  each  pole  of  the  ephemeral  region  must  have  the  same  flux  as  a 
network  element;  (2)  the  number  density  of  network  elements  was  taken 

to  be  independent  of  the  average  f ieldstrength  (2  10  km  ) .  Since 

'3 

each  pole  of  a  typical  ephemeral  region  has  about  1.5  10  Mx  (Golub 
et  al.  1977),  Marsh's  estimate  of  D  applies  only  if  the  average 
f ieldstrength  in  the  network  is  about  3  Gauss.  In  fact  the  process 
can  no  longer  be  described  as  a  diffusion  process,  since  the 
"effective  diffusion  constant"  depends  on  the  average  f ieldstrength  of 
the  "diffused"  flux.  It  is  not  clear  how  Eq.(l)  should  be  modified  to 


include  the 

effects 

of  ephemeral 

regions. 

However , 

the 

ephemeral- region 

effects 

are  not  important 

for  the 

rapid  decay 

of 

large  active  regions,  since  the  f ieldstrengths  are  larger  than  3  Gauss 
in  this  case. 


c)  Small-scale  characteristics 


Observations  of  the  merging  and  disappearance  of  opposite 
polarity  fields  in  a  relatively  quiet  region  on  the  sun  have  recently 
been  described  by  Martin  (1983)  .  In  the  quiet  sun,  network  elements, 
and  poles  of  bipolar  ephemeral  regions,  are  seen  to  move  about  more  or 
less  randomly.  Occasionally,  opposite  polarity  elements  collide,  and 
the  flux  in  each  element  is  seen  to  decrease,  until  the  smaller 
element  completely  disappears.  The  question  arises  whether  the  same 
process  could  be  responsible  for  the  disappearance  of  flux  in  large 
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active  regions,  i.e. ,  does  the  flux  in  active  regions  disappear  due  to 
the  random  coalescence  of  small-scale,  opposite  polarity  fields? 

Wallenhorst  and  Howard  (1982)  studied  changes  in  magnetic  flux 
around  the  time  of  sunspot  disappearance.  Since  the  flux  was  measured 
using  a  raagnetograph  which  saturates  at  a  f ieldstrength  of  about  300 
Gauss,  the  sunspot  fields  are  severely  underestimated.  Therefore  one 
should  expect  to  see  an  increase  in  the  apparent  flux  of  the  active 
region  as  the  sunspot  initially  dissolves.  However,  the  opposite  was 
found:  there  is  a  decrease  of  flux  at  or  near  the  time  the  sunspot 
vanishes.  This  indicates  that  the  flux  lost  by  the  sunspot  disappears 
from  the  solar  surface  on  a  timescale  of  only  1  or  2  days. 

Wallenhorst  and  Topka  (1982)  pointed  out  that  the  disappearance 
of  flux  does  not  seem  to  take  place  at  the  neutral  line  of  the  active 
region.  Rather,  the  magnetic  features  that  move  away  from  the 
dissolving  sunspot  travel  to  the  surrounding  network  across  a 
relatively  field-free  "moat"  (cf.  Sheeley  1969;  Harvey  and  Harvey 
1973) .  When  they  reach  the  outer  boundary  of  the  moat  the  motion 
stops  and  the  flux  elements  seem  to  gradually  fade  away  in  place. 
There  is  no  evidence  for  motion  of  these  features  toward  the  neutral 
line  of  the  region,  which  is  far  removed  from  the  part  of  the  moat 
boundary  where  the  flux  elements  disappear.  The  flux  released  in  the 
break-up  of  the  spot  apparently  does  not  end  up  in  the  surrounding 
network. 

These  observations  suggest  that  the  disappearance  of  flux  in 
large  active  regions  is  not  due  to  the  random  intermixing  of  opposite 
polarity  fields:  flux  disappears  before  intermixing  can  take  place, 
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consistent  with  our  analysis  of  the  diffusion  model  (see  above) .  In 
the  next  section  a  possible  explanation  of  this  phenomenon  will  be 
discussed. 


d)  Summary 


We  suggest  that  there  are  two  different  processes  by  which  flux 
disappearance  takes  place.  One  process  occurs  over  long  periods^time, 
and  is  relevant  to  the  later  stages  of  dispersal  and  decay  of  large 
active  regions.  It  is  caused  by  the  intermixing  of  opposite  polarity 
fields  induced  by  supergranulation  and  granulation.  On  length-scales 
large  compared  with  the  size  of  supergranules,  this  process  can  be 
reasonably  well  described  as  diffusion  (Leighton  1964) .  On  small 
length-scale  the  flux  disappears  by  submergence.  For  average 
f ieldstrengths  less  than  a  few  Gauss,  ephemeral  regions  may  cause  an 
enhancement  of  the  diffusion  rate. 


However,  another  process  of  flux  disappearance  seems  to  take 
place  as  well.  In  large  active  regions,  part  of  the  flux  is  removed 
by  a  process  that  seems  to  involve  the  annihilation  of  fields  at  a 
distance:  bodily  motion  of  magnetic  features  towards  the  neutral  line 
is  not  observed,  and  flux  seems  to  disappear  "in  place"  (Wallenhorst 
and  Topka  1982)  .  As  a  result  of  this  process,  active  regions  decay 
more  rapidly  than  can  be  explained  with  diffusion  models. 
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3.  FLUX  DISAPPEARANCE 

a)  Submergence 

It  has  been  suggested  (e.g.  Zwaan  1978;  Parker  1984a;  van 
Ballegooijen  1982)  that  flux  disappearance  can  occur  by  submergence  of 
magnetic  field.  When  two  flux  tubes  of  opposite  polarity  move  closer 
together  in  the  photosphere,  the  field  above  the  surface  may 
reconnect,  and  form  a  flux  loop.  As  sketched  in  Figure  1,  this  loop 
may  be  pulled  down  below  the  surface  by  the  tension  in  the  magnetic 
field.  Tension  forces  can  overcome  the  buoyancy  forces  only  if  the 
curvature  at  the  top  of  the  loop  is  sufficiently  large.  Presumably, 
convective  flows  below  the  surface  play  an  important  role  in  bringing 
the  flux  elements  closer  together,  thereby  increasing  the  loop 
curvature. 

As  seen  with  a  high-resolution  vector  magnetograph,  the  process 
of  submergence  would  look  as  follows.  Two  patches  of  opposite 
polarity  field  would  be  seen  to  come  closer  together,  until  they 
finally  merge  and  disappear.  Just  before  the  moment  of  disappearance 
the  magnetic  field  would  be  seen  to  turn  horizontal.  We  suggest  that 
the  disappearance  of  magnetic  elements  in  the  quiet  sun,  observed  by 
Martin  (1983) ,  may  indeed  be  due  to  reconnection  and  submergence.  The 
flux-decrease  observed  by  Martin  begins  when  the  separation  of  the 
magnetic  features  is  on  the  order  of  the  spatial  resolution  (about 
3000  km) .  We  suggest,  however,  that  the  actual  submergence  process 
takes  place  on  an  unresolved  scale.  The  submergence  of  a  magnetic 
loop  through  the  photosphere  requires  that  the  magnetic  curvature 
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1 

force  at  the  top  of  the  loop,  B  /(4TIr),  be  larger  than  the  buoyancy 
a. 

force,  B  /(8uH),  i.e.,  the  curvature  radius  R  must  be  smaller  than 

twice  the  photospheric  scaleheight  H.  Since  the  separation  of  the 

footpoints  at  the  base  of  the  photosphere  is  about  2  R,  the  flux 

elements  must  approach  each  other  to  within  a  distance  of  4  H,  or 

about  500  km,  before  submergence  can  take  place.  Thus,  the  observed 

flux-decrease  at  3000  km  separation  must  be  an  effect  of  the  spatial 

resolution  of  the  instrument,  since  at  this  time  the  buoyancy  is  still 

too  large  to  pull  the  flux  loop  down  through  the  photosphere.  The 

loop  must  also  be  fragmented  before  it  can  disappear:  because  of  the 

minimum  separation  of  500  km,  the  thickness  of  each  fragment  can  be  a 

few  hundred  kilometers  at  most.  Therefore  we  suggest  that  only  small 

18 

amounts  of  flux,  of  order  10  Mx  or  less,  are  involved  in  each  event 
(assuming  an  actual  f ieldstrength  of  1000  Gauss  in  the  photosphere) . 

b)  Emergence  of  reverse- polarity  fields 

One  aspect  of  the  submergence  model  is  that  flux  elements  have  to 
move  towards  each  other  for  flux  disappearance  to  take  place?  the 
process  can  occur  only  near  neutral  lines,  and  in  regions  of  mixed 
polarity.  However,  the  observations  discussed  in  Sect.  2  seem  to 
indicate  that  during  the  decay-phase  of  active  regions  part  of  the 
magnetic  flux  disappears  locally,  without  any  motion  of  the  field 
towards  neutral  lines.  To  explain  these  observations  we  propose  that 
a  different  process  of  flux  disappearance  takes  place  as  well.  This 
process,  sketched  in  Figure  2,  involves  the  emergence  of  a  train  of 
small,  reverse-polarity  bipoles  in  between  the  two  patches  of  opposite 
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polarity  of  the  active  region.  The  separation  of  elements  in  each 
bipole  increases  with  time,  bringing  opposite- polarity  elements  of 
neighboring  bipoles  closer  together.  The  magnetic  field  above  the 
surface  then  adjusts  itself  by  opening  up  the  old  bipoles,  and  forming 
new  bipoles  that  connect  the  approaching  elements.  The  first  and  last 
bipole  in  the  train  are  forced  to  reconnect  with  the  original  patches 
of  field.  Subsequently,  the  newly  formed  bipoles  disappear  by 
submergence;  in  the  process  a  fraction  of  the  original  field  also 
disappears.  It  is  important  to  note  that,  in  the  process  described 
above,  there  is  no  motion  of  opposite  polarity  elements  over  large 
distances  across  the  solar  surface.  The  combined  effect  of  the 
motions  within  all  the  bipoles  (emergence,  reconnection  and 
submergence)  is  the  same  as  the  '-ransport  of  flux  over  a  large 
distance,  but  a  single  magnetic  element  does  not  travel  over  this 
large  distance.  As  a  result,  flux  seems  to  disappear  in  place. 

The  emergence  of  reverse-polarity  bipoles  in  decaying  active 

regions  has  not  been  observed.  Thus,  if  the  process  indeed  occurs,  it 

must  take  place  on  a  length-scale  at  or  below  the  spatial  resolution 

of  the  instruments;  the  sensitivity  of  the  magnetographs  also  plays  a 

role.  These  observational  constraints  limit  the  possible  size  of  the 

bipoles  to  about  1000  km,  and  the  maximum  allowable  f ieldstrength  on 

this  scale  is  about  10  gauss.  The  following  calculation  shows, 

however,  that  a  sufficient  amount  of  flux  could  emerge  in  such  bipoles 

to  acount  for  the  rapid  decay  of  active  regions.  A  10  gauss  field  on 

17 

a  scale  of  1000  km  corresponds  to  a  flux  of  10  Mx  per  bipole.  We 
assume  that  every  granule  in  the  active  region  carries  a  new  bipole  to 
the  surface,  and  that  each  bipole  is  part  of  a  train  of  bipoles, 
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connecting  the  two  opposite  polarities  in  the  region  (cf.  Figure  2). 
Tne  width  of  an  active  region,  in  the  direction  perpendicular  to  the 
line  connecting  opposite  polarities,  is  typically  50,000  km,  and  since 
the  width  of  each  bipole  is  about  1000  km,  there  are  of  order  50 
trains  present  at  each  time.  The  granular  turnover  time  is  about  10 

17 

minutes,  so  that  the  bipoles  can  remove  flux  at  a  rate  of  5  10  ' 

2o 

Mx/min,  or  7  10  Mx/day.  Thi&  is  comparable  with  the  decay  rate  of 

to 

5  10  Mx/day  found  by  Wallenhorst  and  Topka  (1982) . 

c)  A  model  involving  sub-surface  reconnection 


To  show  that  our  proposed  model  of  flux  disappearance  in  not 
entirely  unreasonable  we  now  consider  the  possible  origin  of  the 
bipoles.  We  suggest  they  could  result  quite  naturally,  if  the  two 
opposite  polarities  of  the  active  region  reconnect  below  the  surface, 
forming  a  U-shaped  magnetic  flux  tube  in  the  convective  zone  (see 
Figure  3) .  When  the  lower,  horizontal  part  of  this  U-loop  floats  to 
the  surface,  it  emerges  as  a  series  of  bipoles,  that  have  an 
orientation  opposite  to  that  of  the  active  region  as  a  whole.  The 
interaction  of  the  bipoles  with  the  overlying  horizontal  field  of  the 
active  region  leads  to  a  second  phase  of  reconnection,  occuring  in  the 
chromosphere  or  low  corona.  This  reconnection  process  creates  small 
closed-loop  structures  (Figure  3d),  that  may  be  pulled  below  the  solar 
surface  by  granules,  and  may  ultimately  be  destroyed  inside  the 
convective  zone.  In  the  following  we  describe  this  model  in  more 
detail . 
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We  suggest  that  giant-cell  convection,  occuring  in  the  deep 
layers  of  the  solar  convective  zone  (Simon  and  Weiss  1968) ,  is 
responsible  for  driving  part  of  the  opposite- polarity  fields  of  the 
active  region  closer  together,  thereby  forcing  the  fieldlines  to 
reconnect  (cf.  Figure  3b).  This  process  could  take  place  in  a  number 
of  ways.  One  possibility  is  that  part  of  the  active- region  flux 
remains  in  the  convective  updraft  that  was  originally  responsible  for 
the  formation  of  the  active  region.  The  converging  horizontal  flows 
at  the  base  of  the  updraft  push  the  vertical  legs  of  flux  loops  toward 
the  center  of  the  updraft  region.  There  turbulent  or  twisting  motions 
lead  to  the  crossing  of  flux  tubes,  resulting  in  reconnection  of 
opposite  polarity  fields.  This  process  takes  place  on  the  typical 
timescale  of  giant-cell  convection  (about  10  to  30  days) .  Since  the 
giant  cells  apparently  do  not  penetrate  to  the  surface  (Howard  and 
LaBonte  1980;  LaBonte,  Howard,  and  Gilman  1981)  the  flux  tubes  are 
pushed  together  only  in  the  deeper  layers  of  the  convective  zone, 
leading  to  the  formation  of  a  U-loop. 

Once  a  U-loop  has  formed,  the  effects  of  magnetic  tension  and 
buoyancy  act  together  in  pulling  the  loop  upward.  However,  the 
density  in  the  lower,  horizontal  part  of  the  loop  is  high,  and  so  it 
contains  a  large  amount  of  mass.  This  mass  prevents  the  U-loop  from 
freely  escaping  through  the  solar  surface,  since  gravity  will  pull  the 
horizontal  part  back  down  as  soon  as  it  emerges  above  the  photosphere. 
In  order  for  the  U-loop  to  escape  it  has  to  loose  its  mass.  This  is 
difficult  in  view  of  the  high  conductivity  of  the  solar  plasma:  the 
magnetic  field  is  nearly  "frozen"  into  the  gas.  In  fact,  as  long  as 
the  flux  loop  is  slender,  Ohmic  diffusion  can  only  cause  a  further 
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increase  of  the  mass  in  the  loop,  since  pressure  equilibrium  requires 
that  the  gas  pressure  inside  the  loop  is  lower  than  the  external 
pressure  by  an  amount  B  /8X  .  Thus,  Ohmic  diffusion  can  only  drive 
more  mass  into  the  flux  loop.  It  seems  clear  that  U-shaped  magnetic 
loops,  floating  in  the  convective  zone,  cannot  simply  move  out  of  the 
sun. 

To  overcome  this  difficulty  we  propose  the  following  scenario  for 
the  evolution  of  a  U-loop.  As  the  bottom  part  of  the  U-loop  rises 
upward,  the  density  in  this  horizontal  part  decreases  dramatically. 
Since  the  length  of  the  loop,  and  the  mass  inside  it,  remain 
approximately  constant,  the  cross-sectional  area  increases.  This 
results  in  a  decrease  of  the  f ieldstrength,  proportional  with  the 
density.  We  assume  that  the  reconnection  occurs  at  a  depth  of  about 
150,000  km  below  the  surface,  50,000  km  above  the  base  of  the 
convective  zone.  When  the  U-loop  rises  to  a  depth  of  20,000  km  below 
the  surface,  just  below  the  layer  of  supergranules,  the  density,  and 
therefore  the  f ieldstrength,  decrease  by  a  factor  40  (see  the 
convection-zone  model  of  Spruit  1976)  .  If  the  f ieldstrength  was  of 
order  10  Gauss  at  the  reconnection  point,  then  it  is  only  250  Gauss 
at  20,000  km,  small  compared  with  the  local  equipartition 
f ieldstrength  (B^  =  4400  Gauss).  At  this  point  the  supergranular 
flow  starts  to  dominate  the  motion  of  the  field.  We  suggest  that  the 
supergranules  pull  stitches  of  weak  field  up  toward  the  surface;  in 

Figure  4a  we  have  sketched  the  topology  of  the  field.  Each  cell  may 

10 

pull  a  flux  tube  of  typical  size  10,000  km  (flux  of  2.5  10  Mx)  out 
of  the  underlying  horizontal  field.  This  flux  does  not  immediately 
appear  at  the  solar  surface  since  the  f ieldstrength  is  so  low:  the 
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convective  flow,  and  not  buoyancy,  dominates  the  motion.  A  small 
fraction  of  the  flux  reaches  a  depth  of  about  1000  km  below  the 
surface.  Since  the  flux  tubes  are  strongly  distorted  by 
super granulation,  gas  can  flow  down  to  the  lowest  parts  of  the 
fieldlines,  and  so  most  of  the  mass  remains  at  a  depth  of  20,000  km. 
Because  of  these  downflows,  the  f ieldstrength  at  1000  km  depth  is 
difficult  to  estimate;  we  suggest  it  may  be  on  the  order  of  10  to  100 
Gauss.  (Note,  however,  that  the  f ieldstrength  would  have  been  much 
smaller,  B  ;£  0.2  Gauss,  if  the  field  were  to  remain  horizontal  all 

the  way  up  to  a  depth  of  1000  km.) 

Then  a  similar  process  takes  place  due  to  granulation,  i.e.,  on  a 
length-scale  of  1000  km.  As  sketched  in  Figure  4b,  granulation  pulls 
small  flux  tubes  out  of  the  field  below,  and  these  tubes  emerge  at  the 

surface  in  the  form  of  small  bipoles.  The  flux  in  each  bipole  is  a 

17 

few  times  10  '  Mx.  The  bipoles  interact  with  the  overlying  coronal 

field  by  a  reconnection  process  that  occurs  at  X-type  neutral  points 

in  the  chromosphere  or  low  corona  (indicated  by  crosses  in  Figure  4b). 

✓ 

The  reconnection  proceeds  very  rapidly  because  of  the  high  Alfven 
velocity  in  the  chromosphere  and  corona.  The  reconnection  produces 
closed  flux  loops  (see  Figure  4c).  The  smaller  loops  have  a  typical 
size  of  about  1000  km,  but  some  larger  ones,  that  are  rooted  at  about 
20,000  km  depth  and  contain  most  of  the  mass,  also  exist.  We  suggest 
that  the  closed  loops  are  subsequently  pulled  below  the  surface  by  the 
granules.  The  timescale  of  emergence,  reconnection,  and  submergence 
is  on  the  order  of  the  granular  turnover  time,  about  10  minutes.  The 
net  effect  of  the  process  is  the  disappearance  of  field  from  the 
surface. 
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The  rate  at  which  the  horizontal  fields  above  and  below  the 
surface  can  annihilate  each  other  is  determined  by  the  speed  at  which 
granulation  and  supergranulation  can  transport  field  to  the 
photosphere.  A  simple  model  of  this  process  is  that  the  convection 
acts  as  a  diffusive  medium,  which  transports  a  "mean"  horizontal  field 
in  the  vertical  direction.  The  timescale  X  of  the  diffusion  process 

is  determined  by  the  diffusion  constant  *7  ,  and  by  the  thickness  d_  of 

x 

the  convecting  layer:  X  -  /Tj  .  Taking  c[  =  20,000  km,  and  assuming 
that  1 1  is  equal  to  the  diffusivity  for  horizontal  motions  across  the 
surface  (D  =  200  km  /s,  see  Sect.  2) ,  we  obtain  X  m  23  days.  This 
crude  estimate  suggests  that  U-loops  can  disappear  on  a  timescale  of 
one  month. 

To  summarize,  we  have  suggested  that  part  of  the  flux  in  large 
active  regions  is  reconnected  below  the  surface  as  a  result  of  the 
motions  within  giant  cells.  The  U-loops  that  are  formed  in  this 
process  float  to  the  surface,  and  erupt  there  in  the  form  of  small 
bipoles.  The  reconnection  of  these  bipoles  with  the  overlying 
horizontal  field  in  the  corona  leads  to  the  break-up  of  the  loops  into 
many  closed-loop  structures.  These  closed  loops  contain  all  the  mass 
of  the  U-loop,  but  do  not  represent  any  net  flux,  as  seen  with 
low-resolution  magnetographs.  We  conclude  that  reconnection,  first  in 
the  deep  layers  and  then  again  at  the  surface,  may  account  for  the 
apparently  local  disappearance  of  flux. 

The  final  disengagement  of  mass  and  magnetic  field  takes  place 
only  by  the  destruction  of  the  closed  loops.  This  process  may  involve 
the  contraction  of  the  loops  onto  themselves,  leading  to  the 

dissipation  of  magnetic  field  at  the  O-type  neutral  point  at  the 
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center  of  each  loop.  The  dissipation  does  not  necessarily  take  place 
within  the  photosphere:  because  of  their  small  size  and  weak  magnetic 
field,  the  closed  loops  may  be  pulled  below  the  solar  surface  by  the 
downdrafts  at  granular  boundaries.  We  suggest  that  the  closed  loops 
contribute  to  a  weak,  turbulent  magnetic  field  in  the  convective  zone, 
that  does  not  represent  any  "mean"  field.  This  turbulent  field  has  a 
length  scale  on  the  order  of  the  granulation,  and  may  have  a  typical 
f ieldstrength  of  a  few  Gauss.  In  the  course  of  time  the  field  is 
dispersed  over  the  sun  by  convective  motions.  The  inner-network 
fields  observed  by  Livingston  and  Harvey  (1975)  may  be  a  manifestation 
of  the  turbulent  field. 
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4.  CONCLUSIONS 


We  have  analyzed  various  observations  of  flux  disappearance,  and 
suggest  that  active  regions  disappear  more  rapidly  than  can  be 
explained  with  diffusion  models.  There  is  evidence  that  the  flux 
disappears  locally,  without  any  apparent  motion  towards  neutral  lines. 
Such  a  process  cannot  be  readily  explained  in  terms  of  the  dynamics  of 
surface  magnetic  fields;  it  is  necessary  to  invoke  rather  peculiar 
processes  occuring  in  the  convective  zone  below.  In  Sect.  3  we  have 
outlined  one  particular  model,  involving  the  reconnection  of  fields  in 
the  deep  layers  of  the  convective  zone.  Although  this  model  is 
speculative,  it  demonstates  that  one  is  forced  to  consider  rather 
strange  processes  to  explain  the  phenomenon  of  "local"  flux 
disappearance. 

What  are  the  observational  consequences  of  our  model?  Of  course, 
high-resolution  magnetograph  observations  should  show  the  presence  of 
small,  reverse-polarity  bipoles  in  active  regions.  To  detect  these 
small-scale  fields,  long  time-sequences  of  magnetograms  with 
sub-arcsecond  resolution  are  required,  and  the  motions  of  magnetic 
features  on  timescales  of  minutes  to  days  should  be  followed. 
Hopefully,  the  Solar  Optical  Universal  Polarimeter  (SOUP)  on  Spacelab 
2  and  Sunlab  will  provide  insights  into  the  process.  The  Coordinated 
Instrument  Package  (CIP)  on  the  Solar  Optical  Telescope  will  be  able 
to  provide  us  with  further  data  of  this  type. 

Another  way  in  which  the  process  may  manifest  itself  is  by  the 

heating  of  the  atmosphere  that  results  from  the  reconnection  of  the 
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bipoles  with  the  overlying  coronal  field.  Unfortunately,  it  is 
unclear  how  much  energy  may  be  involved  in  this  process,  since  only  a 
small  fraction  of  the  magnetic  energy  is  dissipated  in  the  atmosphere. 
Once  the  field  of  the  U-loops  is  reconnected,  granulation  may  carry 
the  field  down  below  the  surface,  where  most  of  the  magnetic  energy  is 
dissipated. 

Finally  we  note  that,  in  the  present  paper,  we  have  specifically 
considered  those  processes  that  can  affect  a  single,  isolated  active 
region,  and  we  have  ignored  interactions  between  active  regions.  The 
reason  is  that  enhanced  flux-decay  seems  to  be  a  general  phenomena 
that  takes  place  also  in  isolated  regions,  i.e.,  the  interaction  with 
surrounding  fields  does  not  seem  to  be  essential.  Interactions 
between  active  regions  may  be  important,  however,  in  the  context  of 
dynamo  theory.  Parker  (1984b),  for  example,  considered  active  regions 
erupting  on  the  same  toroidal  flux  tube,  and  suggested  that  coronal 
re  onnection  of  flux  from  adjacent  active  regions  may  lead  to  the  loss 
of  toroidal  field  from  the  convective  zone.  Apart  from  the  scale  of 
the  phenomena,  Parker's  mechanism  is  similar  to  ours:  the  sections  of 
toroidal  field  in  between  active  regions  are  like  our  U-loops,  and 
they  may  disappear  in  a  similar  way.  However,  in  contrast  with 
Parker's  model,  our  U-loops  are  formed  within  active  regions,  and 
their  disappearance  does  not  lead  to  a  loss  of  toroidal  flux. 
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Vertical  cross-section  of  convective  zone,  illustrating  the 
disappearance  of  flux  by  submergence:  (a)  two  opposite 
polarity  flux  tubes  move  closer  together,  and  (b)  are  pulled 
below  the  surface. 


The  field  distribution  in  the  photosphere  during  the  process 
of  flux  disappearance:  (a)  a  train  of  small  bipoles  erupts  in 
between  two  patches  of  opposite  polarity;  (b)  positive 
(negative)  flux  elements  move  a  small  distance  to  the  left 
(right),  effectively  transporting  flux  over  a  much  larger 
distance. 


Illustrating  the  disappearance  of  flux  by  U-loop  formation: 
(a)  flux  tubes  move  closer  together  in  the  deep  layers  of  the 
convective  zone;  (b)  reconnection  leads  to  formation  of 
U-loop;  (c)  horizontal  part  of  u-loop  floats  up  to  the 
surface,  and  erupts  as  a  series  of  bipoles;  (d)  reconnection 
of  bipoles  with  coronal  field  produces  closed  flux  loops, 
which  then  sink  below  the  surface. 


The  emergence  of  U-loops:  (a)  supergranulation  distorts  the 
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field,  and  mass  flows  to  the  lowest  parts  of  the  loop;  (b) 
granulation  brings  stitches  of  field  above  the  surface,  where 
it  reconnects  with  the  overlying  coronal  field  (crosses  mark 
X-type  neutral  points);  (c)  as  a  result  of  reconnection  closed 
loops  are  formed,  that  are  pulled  below  the  surface  while  new 
flux  emerges. 


b  )i 


